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A rewiew of the development of formalism to describe correlations and fluctuations 
in multiparticle production is presented. The milestones of history of this develop- 
ment and some obvious sources of correlations in experimental data are recalled. 
Some recent subjects concerning the correlations are shortly discussed. 

1 Introduction 

This talk|^| is supposed to give an introduction for the session devoted to 
theory of correlations and fluctuations. Thus the standard form should be a 
review of methods employed recently in the investigations of these phenomena. 
However, such a review was already presented by Hans Eggers in the previous 
sessionn and I would not dare to try to add anything to his comprehensive 
presentation. 

Therefore I am going to review instead the methods and facts which are 
with us for such a long time that most of them are often forgotten. In particu- 
lar, I would like to remind you that this is not only the jubilee 30th Symposium 
on Multiparticle Dynamics (called in the seventies by a few other names), but 
also the 30th anniversary of the beginnings of modern history of investigations 
on correlations and fluctuations in multiparticle production. Thus I shall re- 
call some milestones of this history (obviously this short note is not a serious 
historical review, but rather a highly biased personal account of some facts) . 
Then I review briefly obvious sources of correlations in experimental data and 
conclude with a short discussion of some recent developments. 

2 Highlights from the History of Correlations in Multiparticle 
Production 

In the fall of 1970 Mueller adapted for the processes of multiparticle pro- 
duction a mathematical formalism known already well in the statistics and 
some branches of physicso. The formulae are most transparent for inclusive 
distributions, i.e. for the densities of particles in momentum space averaged 
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over all events. One defines the relations between the one- and many-particle 
densities p q and correlation functions c q by 



C2(P1,P2) = P2{pi,P2) - pi{pi)pi(p 2 ), (1) 

Ca(pi,P2,P3) = P(Pl,P2,P3) - Pl(Pl)c2(P2,P3) - Pi {P2)c 2 (pi , P3) 

~ Pl(P3)c2(Pl,P2) - P{P1)P{P2)P{P3) (2) 

etc., where the densities p q are normalized to the factorial moments of multi- 
plicity distributions 

Pqipu ■■■Pq)d J 'Pi.--d i p q =<n(n -l)...(n- q+ 1) >=< n q > . (■'>) 

Thus 

h = J C2{pi,P2)d 3 p 1 d 3 p 2 =< n 2 > - <n > 2 = D 2 - <n>, (4) 

fa = J c 3 (p 1 ,p2,P3)d 3 pid 3 p 2 d 3 p 3 =< n 3 > -3 < n > f 2 - < n > 3 

= < n 3 > -3 < n >< n 2 > +2 < n > 3 , (5) 

etc. These formulae justify the standard connection between "correlations" 
and "fluctuations" terms: we see that non-Poissonian fluctuations in multi- 
plicity distributions (i.e., /„ ^ for q =/= 1) imply non-zero correlations. In 
a more recent notation, Mueller's "correlation integrals" f q are called "fac- 
torial cumulants" and denoted by K q . It is often convenient to use "scaled" 
quantities divided by a proper power of < n >, i.e. 

F q =< n q > I < n > q , K q = n q / < n > g . (6) 

One distinguishes usually between the "short-range" correlations (SRC), i.e. 
terms in c q which vanish for large differences between pi and integrate in n q -s 
to terms growing with < n > slower than < n > 9 , and the "long-range" cor- 
relations (LRC) yielding K q -s growing as < n > q . Details of such a distinction 
depend usually on the particular model in which they are calculated. One 
of the first profits from the efforts to identify the souucas of LRC was the 
separation of inelastic diffraction as the source of LRCB'Q which dominated 
the early data. 

Relations between the factorial moments < n q > and factorial cumulants 
K q for any values of q may be derived from their expressions by the generating 
function of the multiplicity distribution P(n) 

G(z)^^P(n)z" (7) 
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which are 



< n q >=d q G(z)/dz q \ z=1 , K q = d q lnG(z)/dz q \ z=1 . (8) 

This yields for the scaled moments and cumulants 

Fi = Ki = l, F 2 = l + K 2 , F 3 = l + 3K 2 + K 3l (9) 

etc. Recently it was noted that the ratios of cumulants to moments H q = 
K q /F q are even more useful for testing the models, as the errors largely cancel 
in such ratios and even for q up to 10 the data allow to calculate H q with 
relatively small uncertainties. 

The apparent sharp contrast between the "differential" quantities, as the 
densities p q and correlation functions c q and the "integral" quantities as the 
factorial moments < n q > and cumulants n q is in fact slightly artificial. We 
never measure really "differential" distributions, as all the experimental his- 
tograms are integrals over bins. On the other side, both the experimental 
limitations and theoretical suggestions led to the use of integrals restricted to 
some part of phase-space in the LHS of Eq.(|J) and the following formulae. 

This line of investigations included the subject of forward-backward cor- 
relations, the approximate KNO scalingQ defined by 

I a 

P(n) = /( ) (10) 

v ' <n> J K <n>' y ' 

not only in the whole pase-space, but also in its parts defined by rapidity bins, 
and in particular the dependence of normalized moments on energy and on 
the size of rapidity bin. The negative binomial distribution. (NBD) reinvented 
for multiparticle production by Giovannini and Van Hovecl 

was found to describe reasonably well all the measured multiplicity distribu- 
tions. For a perfect fit n equals the experimental average multiplicity < n >. 
Dependence of the parameters n and k on the energy and bin width allowed 
for important tests of various models. The valuable feature of NBD is that the 
superposition of N NBD-s with parameters n, A; is a NBD with parameters 
Nn, Nk, as one may easily check using the generating function 



G N . NBD (z) = [G n / BD (z)] N = [1 + ^(1- z )]- kN = G N N n d ^ k {z). (12) 



This shows in a particularly obvious way that superposition of N independent 
sources yields a multiplicity distribution with correlations reduced by a factor 
1/N compared to the individual distribution, since for NBD n 2 = 1/fe. 
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The next important step was the proposal to investigate scaled factorial 
moments in the limit of very small rapidity bins 13 (averaged over the bin 
position) 

< f > s Sgi < a >• — . ( 13 ) 

This proposal, originally motivated by a search for possible " intermittency" 
effect due to the phase transition, resulted indeed in the discovery of expected 
power-like increase of < F q > with M (or the inverse bin width). However, 
the very universality of the effect makes its interpretation rather difficult. 
In any case, it allows to investigate in detail many interesting short-range 
correlation effects absent in the standard Monte Carlo generators, and in par- 
ticular the Bose-Einstein (BE) interference effects, which seem to dominate 
the observed increase of moments for hadronic collisions. Later, 2-D anQ 
D bins in momentum space were used instead of simple rapidit^bins, 
and variables guaranteeing uniform distributions were introducedHtil. Most 
important modification, however, was the replacement of the sum over "rect- 
angular" bins by the " snake" integrals of densities! 2 ! over phase-space regions 
defined by small jnomentum differences, which suppressed greatly the statis- 
tical fluctuations!^. In particular, the ratios defined for investigating the BE 
effects were rediscovered 

fl( g 2 j = II d 3 P id 3 P2 p 2 ( Pl ,p 2 )S(Q 2 + (p! - P2 ) 2 ) 
J J d 3 pid 3 p 2 p 1 (p 1 )p 1 (p 2 )S(Q 2 + (pi -P2) 2 ) 

and many interesting ideas on possible emergence of intermittency in such 
ratios and other quantities were presented. Many of these and other results 
may be found in more detail in a review by De Wolf, Dremin and Kittcltil. 



3 Some Obvious Sources of Correlations 

There are many well known sources of correlations in multiple production. We 
recall some of them here to remind the reader that any analysis of data which 
neglects them may be very misleading. In particular, the model parameters 
fitted to the data may change significantly if these effects are " switched off" . 
Two such "obvious" effects are indisputable: 

i) resonance effects. It is known that the majority of observed pions comes 
from resonance decays. This induces strong positive short-range correlations 
for unlike-sign pions and weaker but non-negligible correlations for like-sign 
pions. 
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ii) energy-momentum conservation. It gives negative correlations which ex- 
tend over full phase-space, but integrate to terms in k^-s which grow slower 
than < n > q (thus they are neither typical SRC, nor LRC). 
Other " obvious" sources are model-dependent but occur for such large classes 
of models that many people regard them as evident. Let us give two examples: 
i) string effects. They are common for all the models which assume that mul- 
tiple production is the result of colour separation and the subsequent string 
breaking. Among others, they include effects of flavour ordering, local com- 
pensation of transverse momenta, long-range effects from flavours at the string 
ends and the inter-string effects for multi-string events (" drag effect" , deple- 
tion in "no-string" sectors between jets). 

il) BE (or Fermi-Dirac) interference effects. Their very existence is certain, 
but their impact on correlation data depends on extra assumptions. They 
induce positive very short-range correlations for identical pions and other 
mesons, and negative SRC for identical baryons. For the better known BE 
effects one may investigate not just the standard "source size" effects, but 
also spatial asymmetry of the source and higher order effectsEj. 

Obviously, many of these effects cannot be introduced by any analytic 
formulae. Thus it is very useful to build the model ideas into Monte Carlo 
generators, in which all such effects may be incorporated. Otherwise no con- 
clusions about the agreement or disagreement between the models and data 
are fully reliable. 

4 Some Recent Developments 

Most important development of the last decades was the possibility to in- 
vestigate the correlations in the framework of QCD, and not just in more or 
less ephemeral models. This is described in widely accessible reviewalil and 
we will not recall it here. However, let us mention that recently the ideas 
of using infrared-safe quantities^ seem almost forgotten (perhaps as insuf- 
ficiently ambitious?). One tries to derive predictions for the experimentally 
measured quantities (as the factorial moments and cumulants) not just in the 
lowest order of pQCD, but solving numerically the equations for generating 
functionstll. One of the most interesting physical problems seems to be the 
role of QCD interference effects in the correlation data and in particular in 
the multiplicity distributions. These and other subjects are discussed in a 
recent review by Dremin and GaryEa. 

The most lively discussed experimental data concern the BE effects. Wait- 
ing for the flood of RHIC results one develops many methods which are already 
used successfully at LEP. Again, we refer to the recent review by Weinenij 
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for more details. Here let us mention only two subjects which are discussed 
in this symposium: the controversy on the interference effects betaieen pions 
from the decay of two different W-s, which is still not resolved i3 and the 
proposals to treat the effeci-in e + e" collisions in a completely different way 
than in heavy ion collisionsE2L [Xhe possibility to model the details of the effect 
(e.g. the source asymmetryE!rE3) in MC generators is also widely debated. 

Summarizing, there is a lot of work to be done experimentally and theo- 
retically to exploit fully the potential of the correlation phenomena in multi- 
particle production. The aim of this mini-review was to recall some highlights 
from the history of this subject and to suggest that we should not forget the 
old results, which may be as relevant today as they were in seventies. This 
seems to be especially true for the BE effect, where the interpretation of the 
data seems today equally difficult as it was thirty years ago. 
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